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ANGLE-OF- ATTACK-MODULATED TERMINAL POINT 
CONTROL FOR NEPTUNE AEROCAPTURE 

Eric M. Queen* 

An aerocapture guidance algorithm based on a calculus of variations ap- 
proach is developed, using angle of attack as the primary control variable. 

Bank angle is used as a secondary control to alleviate angle of attack ex- 
tremes and to control inclination. The guidance equations are derived in 
detail. The controller has very small onboard computational requirements 
and is robust to atmospheric and aerodynamic dispersions. The algorithm 
is applied to aerocapture at Neptune. Three versions of the controller are 
considered with varying angle of attack authority. The three versions of 
the controller are evaluated using Monte Carlo simulations with expected 
dispersions. 

INTRODUCTION 

“Aerocapture” is the use of aerodynamic forces to capture a satellite into an elliptic orbit from a 
hyperbolic approach orbit. Because aerocapture promises great reductions in mass needed to achieve 
orbit about certain planets, it has been proposed for several advanced missions. [1, 2, 3] Often, the 
vehicles used for these missions are blunt bodies with a small, nearly constant lift-to-drag ratio 
(L/D). This approach has heritage dating back to the Apollo era. [4] The blunt body configurations 
considered have been designed to fly at a constant angle of attack. The trajectories flown by these 
vehicles have been controlled by modulation of the bank angle. Bank angle magnitude is used to 
control the apoapsis altitude and sign of the bank angle is used to control the inclination or wedge 
angle. 

Recently, a mission has been proposed for aerocapture at Neptune. Neptune’s deep gravity well 
leads to very high entry velocities, on the order of 30 km/s. For these high entry velocities, bank 
angle modulation may not provide sufficient control authority' to ensure accurate aerocapture. To 
overcome this issue, a combination of angle of attack modulation and bank angle modulation may 
be required. 

The guidance algorithm for the Apollo Earth return used a terminal point controller based 
on calculus-of-variations theory' for the final phase of flight [5]. Later, this approach was applied 
to aerocapture using bank angle modulation as the control [6, 7]. The current work develops a 
terminal point controller for aerocapture with angle of attack as the control variable. This type of 
controller is desirable for several reasons. The first is small computational requirements. The current 
engineering code is about 350 lines for the onboard portion of the algorithm, with no iterative loops. 
Also, for the types of dispersions and uncertainties expected, experience has shown the algorithm 
to be relatively robust. Perhaps most inportanth , when the algorithm is pushed beyond reasonable 
limits, it degrades gracefully'. 

The General Theory' section gives details of the mathematical formulation of a two-point bound- 
ary' value problem. The section on Application to Aerocapture applies this process to the problem of 
aerocapture. The Mission section gives a brief description of the Neprtune Aerocapture mission. The 
Numerical Results section shows results of numerical Monte Carlo simulations that were performed 
to establish performance of the algorithm. 
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GENERAL THEORY 


The general theory of terminal point control is widely available. See, for instance, Applied Optimal 
Control by Bryson and Ho[5]. The following development follows Bryson and Ho. 

Consider a system described by differential equations 

X = f(x, u, t) (1) 


with x(tO) given. For aerocapture, it is convenient to write the cost function as: 

J = <t>(x(tf),tf) 

after adjoining the differential equations 1 to equation 2 and integrating by parts: 

J = &(x(tf),tf) - X T (tf)x(tf) + A r (io)^(io) 

r t. 


/to 

Next consider the first order variations in J 


J |a t (t)f(x,u,t) + A T (t)x(£) 


dt 


[A t 5x] 


"-[(MH- 

+ £[( aT 5 ; +f ) 5 i+ ( aT ^)H‘ ! ‘ 


A convenient choice for A is: 
with boundary conditions 

which yields 


A =- iT r 


\ T (tr) = 

{/) dx(t f ) 


5J = [A T H t _t 0 + 1' [*%) ««] dt 

Particular solutions to the differential equations are:[6] 

wA xTMh{, ° )+ n(^) 

A (t) = J & dt 


Sx(t ) = 


<5^ dt 


and 


From equations 7 and 8 it follows that 


5J = X T (t)5x(t) 


(2) 


(3) 


(4) 

(5) 

( 6 ) 

(7) 

( 8 ) 

(9) 

( 10 ) 


Which means that we have linearized the change in cost function with respect to deviations in the 
states. The coefficients in the linear equations are determined by equations 5 and 6. 
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APPLICATION TO AEROCAPTURE 

In-Plane control 




The relevant performance index J for aerocapture is AV , the change in velocity required to reach the 
desired orbit. A nominal trajectory will be determined which has a value of AV that is acceptable 
to the mission (preferably small) a guidance which drives 8J to zero will maintain AV near the 
nominal value. 

We take as our state 

V(t) 

*(*)= 7(*) | (11) 

L m 


with state equation 


f(x(t),u(t)) = 


~(i+psin(T)) 


V sin(y) 


and costates 


which leads to: 


m= 


Av(*) 
A 7 (t) 
A n(t) 


de 

chz* 

_2D 1 
V* m 

+ (tT+A — Tvh 2 ) cos (t) 

sin( 7 *) 


( 


J- 

V* 


9 cos ( 7 *) 
V* cos( 7 *) 


D*/m 

h a 

L*<x>s( 4 >) . V* 008 ( 7 *) 

mV~h: (r.+k *) 2 

0 


( 12 ) 


(13) 


(14) 


where the superscript * indicates values taken from a nominal trajectory. Taking AV as our perfor- 
mance index, equation 7 becomes 


* 


5AV = [A T 5*] 




■mi 




(15) 


with 


A T (*/) = 


9AF(t/) 

&r(£/) 


(16) 


Which is a two-point boundary value problem. To solve this problem, we set 5AV = 0 and solve for 
the control 5u. Unfortunate!}, the solution for 5u is impossible unless we assume a form for 5u. For 
convenience, we will assume 5u is a constant and a scalar. In which case we have 



we define 

After some manipulation the time derivative of Xu is determined to be: 


d\ u 

dt 


-A T (t) 


am 

du(t) 


(17) 


(18) 


(19) 
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and 

Equation 17 expands to 


K(tf) = o 

j \h.5h + A ySV + \-y8y 

An 


Taking the first partial of h with respect to y yields 


8y = 


5k 


V'cos( 7 ) 


We assume 


D= -p 0 eq>(-h/h a )V 2 C d A 


and take the partial of D/ m with respect to h yields 

. 8{D/m)mh s 
5h = 

Substituting into equation 21 we have 

-\ k S( p/y™ h - + \ v sv + a. - Lh 


5u = 


7 V cos(y) 




( 20 ) 

( 21 ) 


( 22 ) 

(23) 


(24) 


(25) 


The angle of attack is assumed to respond much more quickly than the bank angle to control 
commands. Hence, we will treat angle of attack as the control (ii = a) and bank angle as a slowly 
varying parameter. We will use this parameter to reduce excursions in the angle of attack. 

The equation for the “control costate' ' is then: 


• _ xT df T _ Ay dD A 7 cos(<£*) dL 
da m da mV* da 


(26) 


It is possible that cos(o) is far from cos(o*), so we also include 

Ay dD A-y dL 

° m da mV* da 

which corresponds to cos(o) = 1, and 

Ay dD A 7 dL 

* m da ^ mV* da 


(27) 


(28) 


which corresponds to cos(o) = —1. These three derivatives are used to produce three different values 
for Ac. The value for A,, used in equation 25 is determined by interpolating between the three values 
based on the current value of o. 

The bank angle is modulated to attempt to drive the angle of attack back to its nominal value. 



j*:(180-MI) 


a > 0 
a < 0 


Thus the control equations are: 


ac = a* + 5a 


and 


<t>c = 4>* + 5<p 


where 5<j> is from equation 29 and 5a is from equation 25; 


(29) 

(30) 

(31) 
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Out-of-Plane control 

The preceding discussion only establishes the magnitude of the bank angle. The sign of the bank 
angle is used to control the final inclination or wedge angle. (Inclination is the angle between an 
orbit and the equator, while wedge angle is the angle between two orbits. The difference being that 
i to achieve a zero wedge angle, one must match both inclination and longitude of ascending node.) 

The sign of the bank angle is reversed whenever the inclination or wedge angle exceeds deadband 
limits. The deadbands are specified as a function of velocity, typically with a funnel shape that 
„ narrows down as the exit velocity is approached. A “latch” logic is included so that once a bank 

reversal is initiated, it must be completed before another reversal can be commanded. The current 
work does not impose tight limits on the inclination. 


Design Trajectory 

It is clear from equations 30 and 31 that the algorithm depends strongly on the design trajectory. 
For the lowest AV\ an ideal trajectory would exit the atmosphere full lift down and just meet the 
desired apoapsis. This would make the periapsis as high as possible. However, any atmospheric 
dispersions or vehicle modeling errors would cause the vehicle to fail to meet the required orbit. For 
increased robustness to atmospheric disturbances, the design trajectory has an angle of attack near 
the midpoint of its range and bank angle near 90 degrees. 

MISSION 

The current guidance scheme is applicable to any aerocapture mission with the capability of angle of 
attack modulation. For demonstration purposes, it was applied to a proposed mission of aerocapture 
at Neptune. One of the primary goals of this mission is a flyby of Neptune’s moon, Triton. To match 
orbits with Triton, the vehicle must have an apoapsis of about 430,000 km. This is a particularly 
difficult aerocapture mission, because the atmospheric exit velocity must be so close the to escape 
velocity. 

The preceding guidance algorithm was coded in Fortran and incorporated into the Program to 
Optimize Simulated Trajectories II (POST II) [8]. The atmosphere model used is NeptuneGram 
Version 0. 

The vehicle chosen is a flat-bottomed ellipsled with nominal L/D of 0.8 at an angle of attack of 
40 degrees. The vehicle has an entry mass of 2450.0 kg and a ballistic coefficient of 257.85 kg/?n 2 . 

The entry states have a nominal entry flight path of 12.5 degrees. For Monte Carlo analysis, the 
states were provided at 60 seconds prior to atmospheric entry. The mean velocity was 28.85 km/s 
with a standard deviation of 0.01 km/s. The flight path angle variability 7 was 0.175 deg 1-sigma. 
The axial force coefficient of the vehicle varied by ± 30% and the normal force coefficient varied by 
± 50%. 


NUMERICAL RESULTS 

Nominal Trajectory 

A nominal altitude-velocity curve is shown in figure 1. At the edge of Neptune’s sensible atmosphere 
(taken as 25764 km radius), the escape velocity is 23 km/s. Figure 1 shows how close the exit 
condition is to escape velocity. The control histories from a nominal trajectory are shown in figure 2. 
The angle of attack is assumed to respond instantly to commands, while the bank angle has a 
simulated controller in the loop that limits bank accelerations to 5 deg/ 5 2 and bank rates to 20 
deg/s. 
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Figure 1: Nominal Altitude-Velocity curve. 

Monte Carlo Analysis 

Three Monte Carlo analyses were performed on three different versions of the a-modulated Terminal 
Point Controller. The three controllers differed in the range of angle-of-attack permitted. All three 
had a nominal angle-of-attack of 40 degrees, and the different versions were allowed excursions away 
from that of plus or minus 10, 5, and 3 degrees. As expected, as the allowable angle of attack range 
is reduced, performance is degraded, both in terms of accuracy of apoapsis as well as AV required. 

Figure 3 shows the AV required to meet the target orbit versus inclination at atmosphere exit 
for 2000 Monte Carlo cases with and angle of attack range of ± 10 degrees. Figures 4 and 5 show 
the same results for angle of attack ranges of 5 and 3 degrees respectively. In Figure 3 the control of 
AV is very good most cases are less than 150 m/s with only a few above 200 m/s. When the angle 
of attack range is reduced to 5 degrees, Figure 4 shows that the AV values increase, some as high 
as 600 m/s. When the angle of attack range is further reduced to 3 degrees, Figure 5 shows that 
many of the AV values exceed 800 m/s. 

Figure 6 shows the maximum sensed acceleration for each of the three monte carlo simulations. 
As can be seen the maximum acceleration increases as the angle of attack range increases. This 
increase is expected due to the higher lift force associated with higher angle of attack. 


CONCLUSIONS 

A guidance algorithm for aerocapture using angle of attack modulation as the primary control has 
been presented. Limitations have been imposed on the range of angle of attack allowed by the 
algorithm. Three different angle of attack ranges have been examined via Monte Carlo analysis. It 
has been shown that as the angle of attack range increases the performance of the guidance algorithm 
improves, but the sensed acceleration also increases with angle of attack range. 
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Figure 2: Nominal Control profiles. 
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Figure 3: Ten degree Angle of Attack Range Monte Carlo. 
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Figure 4: Five degree Angle of Attack Range Monte Carlo. 













